High‐accuracy practical spline‐based 3D and 2D integral transformations in potential‐field geophysics     

Bingzhu Wang  Stephen S. Gao  Kelly H. Liu  Edward S. Krebes 《Geophysical Prospecting》2012,60(5):1001-1016

Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.  相似文献   

Gravity field and tectonic features of Block L2 in the Lamu basin,Kenya     

Bingqiang Yuan  Wensheng Xie  Guihe Liu  Chunguan Zhang 《Geophysical Prospecting》2012,60(1):161-178

A gravity survey on the scale of 1: 250 000 was carried out in Block L2 located in the Lamu basin of south‐east Kenya in order to study tectonic features and find out favourable petroleum prospects in the block. This paper, through data processing and synthetic interpretation of the measured gravity data in the block, discusses characteristics of the gravity field and their geological implications, determines the fault system and the basement depth, analyses features of the main strata, divides structure units and predicts favourable petroleum zones. In the block, the regional gravity anomaly is mainly caused by the inclined Moho surface that rises in the east and subsides in the west topographically and the Bouguer gravity anomaly primarily reflects the superimposition of the gravity effect derived from the Moho surface and the basement relief. Two groups of faults extending NW (NWW) and NE (NEE) respectively are dominant in the block and their activities resulted in the framework of east‐west zoning and south‐north blocking. The basement depth greatly changes in an alternative high and low pattern. The Permian‐Triassic, Jurassic and Tertiary strata are extensively developed, while the Cretaceous is only developed in the east of the block. Structurally, the block can be divided into five units, of which the Tana sag shows excellent source‐reservoir‐seal associations and is a favourable target for future petroleum exploration.  相似文献   

球体位错理论计算程序的总体设计与具体实现   总被引：3，自引：0，他引：3  

付广裕  孙文科 《地震》2012,32(2):73-87

本文详细介绍了球体位错理论计算程序的总体设计思想、 各类配套文件的具体内涵以及各类输出文件的物理含义, 同时介绍了程序的使用方法和注意事项, 为读者独立使用该程序提供参考。 球体位错理论计算程序主要由三部分组成: ① 位错格林函数计算程序, 基于具体的球对称地球模型提供离散的二维格林函数数值框架; ② 积分计算程序, 对离散的格林函数数值框架进行双二次样条插值运算, 并对四类独立位错源对应的格林函数进行适当组合, 从而计算出任意位置任意类型震源在地表产生的同震变形(含位移、 应变、 重力变化和大地水准面变化); ③ 辅助文件, 用于提供发震断层模型和计算点位信息。 一般情况下, 读者不需要理解位错格林函数计算程序和积分计算程序, 只需要对辅助文件提供的信息进行针对性改动, 就可以计算目标地震在目标观测站引起的同震变形。  相似文献   

多网重力数据联合平差和系统误差改正   总被引：1，自引：0，他引：1  

李真  陈石  秦建增  卢红艳  徐伟民 《地震》2012,32(2):95-104

华北地区现行的地震重力监测网由于独立联测, 自成体系, 绝对控制点较少, 获得的重力场动态变化图像在绝对控制较弱的地区会产生畸变, 也不能有效的消除仪器标定系统所引起的测量误差。 本文以华北大网作为控制网, 结合绝对重力点进行联合平差计算, 使各省区地震重力监测网平差精度提高了(0~4)×10^-8 m/s², 平差起算基准相统一, 绝对控制较弱的地区畸变减小。 采用相关系数检验法, 对华北地区各省区地震重力监测网2010年二期重力平差计算结果进行检验, 部分重力点位系统误差高达60×10^-8 m/s², 经过系统误差改正后, 重力场动态变化趋于一致。 经过各项改正后, 华北地区2010年重力场动态差分图像对地震前兆具有更好的反映。  相似文献   

东海陆架盆地雁荡低凸起综合地球物理解释及其成因探讨   总被引：4，自引：0，他引：4       下载免费PDF全文

杨传胜  李刚  栾锡武  杨长清  龚建明  杨艳秋  索艳慧 《地球物理学报》2014,57(9):2981-2992

通过对东海陆架盆地西部地震和重磁资料的综合地球物理解释,对雁荡低凸起展布形态进行了细致刻画,凸起呈NE方向不连续展布于瓯江凹陷和闽江凹陷之间,长约170 km、宽约15~50 km.地震资料揭示雁荡低凸起上广泛发育了侏罗纪与白垩纪地层,厚度约为500~1500 m,展布面积约5000 km²,局部缺失中生界地层.凸起两侧中生代盆地结构差异明显,西侧瓯江凹陷为典型的断陷盆地,东断西超、断裂发育,半地堑、掀斜断块等中生界构造样式发育;东侧闽江凹陷为坳陷型盆地,断裂、火成岩不发育,挤压背斜、断背斜、反转构造等中生界构造样式发育.自由空间重力异常图与剩余重力异常图上,凸起表现为一系列NE向团块状重力高值区,而磁力异常ΔT图上则表现为深部磁场强度低的特征,火成岩影响部位可见磁力高值异常.综合凸起及邻域重磁震、莫霍面深度等地质地球物理资料,认为雁荡低凸起为一元古界组成的古隆起,区别于东部的台北低凸起.同时,结合区域构造演化及沉积特征,推测侏罗纪时期雁荡低凸起与浙闽隆起区可能连为一体,晚白垩世近东西向伸展作用下浙闽隆起发生裂陷进而形成了雁荡低凸起.  相似文献   

基于重力地质法的南中国海海底地形反演   总被引：3，自引：0，他引：3       下载免费PDF全文

欧阳明达  孙中苗  翟振和 《地球物理学报》2014,57(9):2756-2765

根据重力地质法（GGM）,利用南中国海海域内63179个船测控制点水深将测高自由空间重力异常划分为长波参考场和短波残差场,并反演出了该海域112°E-119°E,12°N-20°N范围的1’×1’海底地形模型,该过程中使用的海水和海底洋壳密度差异常数1.32 g·cm^-3通过实测水深估计得到.利用反演得到的GGM模型对剩余的10529个检核点船测水深插值计算后与实测水深进行比较,其较差结果的均值为-1.64 m,标准差为76.95 m,相对精度为4.06%.此外,根据船测点数量、分布和海底地形的不同,选择了三个海域进行统计,结果表明：在船测控制点分布均匀的海域,GGM模型精度优于ETOPO1模型,在控制点过于分散的海域其精度会有所下降,但好于船测水深的直接格网化结果.为进一步探究检核点的较差结果中出现较大数值的成因,本文对精度较差的点位进行了单独分析,选择了两条船测航迹剖面进行了研究,并分析了检核点的水深较差、相对精度与水深和重力异常的关系,结果表明：GGM模型精度受水深和重力异常的相关性影响较小,受海底地形复杂程度影响较大,地形坡度变化平缓海域的预测精度明显高于海山地区.最后,综合GGM模型和ETOPO1模型优势,利用所有船测水深作为控制,生成了综合的海底地形模型.  相似文献   

A theoretical study on terrestrial gravimetric data refinement by Earth gravity models     

Mehdi Eshagh 《Geophysical Prospecting》2014,62(1):158-171

The idea of this paper is to present estimators for combining terrestrial gravity data with Earth gravity models and produce a high‐quality source of the Earth's gravity field data through all wavelengths. To do so, integral and point‐wise estimators are mathematically developed, based on the spectral combination theory, in such a way that they combine terrestrial data with one and/or two Earth gravity models. The integral estimators are developed so that they become biased or unbiased to a priori information. For testing the quality of the estimators, their global mean square errors are generated using an Earth gravity model08 model and one of the recent products of the gravity field and steady‐state ocean circulation explorer mission. Numerical results show that the integral estimators have smaller global root mean square errors than the point‐wise ones but they are not efficient practically. The integral estimator of the biased type is the most suited due to its smallest global root mean square error comparing to the rest of the estimators. Due largely to the omission errors of Earth gravity models the point‐wise estimators are not sensitive to the Earth gravity model commission error; therefore, the use of high‐degree Earth gravity models is very influential for reduction of their root mean square errors. Also it is shown that the use of the ocean circulation explorer Earth gravity model does not significantly reduce the root mean square errors of the presented estimators in the presence of Earth gravity model08. All estimators are applied in the region of Fennoscandia and a cap size of 2° for numerical integration and a maximum degree of 2500 for generation of band‐limited kernels are found suitable for the integral estimators.  相似文献   

北京地区主要活断层活动参数的重力反演     

董运洪  刘峡  韩月萍  杜雪松  塔拉 《地震研究》2014,37(3)

采用1991~1999年首都圈重力测网的重力测量资料,根据地质调查结果确定断层初始运动模型,利用稳健—贝叶斯最小二乘方法反演北京地区4条活断层的活动参数(走滑、倾滑、引张)。反演结果表明:在1991~1999年间,北京地区活断层的运动性质均以走滑和垂直升降为主,兼有伸展特征,但运动速率都不大,约为2 mm/a。北京地区北西走向断层反演的运动特征为左旋走滑正断层,北东走向断层反演的运动特征为右旋走滑正断层。  相似文献   

Measured and modelled absolute gravity changes in Greenland     

《Journal of Geodynamics》2014

In glaciated areas, the Earth is responding to the ongoing changes of the ice sheets, a response known as glacial isostatic adjustment (GIA). GIA can be investigated through observations of gravity change. For the ongoing assessment of the ice sheets mass balance, where satellite data are used, the study of GIA is important since it acts as an error source. GIA consists of three signals as seen by a gravimeter on the surface of the Earth. These signals are investigated in this study. The ICE-5G ice history and recently developed ice models of present day changes are used to model the gravity change in Greenland. The result is compared with the initial measurements of absolute gravity (AG) change at selected Greenland Network (GNET) sites.We find that observations are highly influenced by the direct attraction from the ice and ocean. This is especially evident in the measurements conducted at the GNET station near the Helheim Glacier. The effect of the direct attraction diminishes at sites that are more than one degree from the source. Here, the dominant signal is the effect of the elastic signal from present day ice mass changes. We find agreement between the measured and modelled gravity changes at all but one site. This agreement only holds when the direct attraction is considered. For one site, there is no agreement, indicating that some improvements to the modelling results or the processing of the gravity data are needed. In addition, more AG measurements are needed to strengthen the time series of gravity change.  相似文献   

关于黄海深部构造的地球物理认识   总被引：11，自引：7，他引：4       下载免费PDF全文

郝天珧  黄松  徐亚  李志伟  张丽莉  王家林  Suh Mancheol  Kim Kwangkee 《地球物理学报》2010,53(6):1315-1326

利用黄海海区重磁观测数据,以最新的地震层析成像和浅层反射地震探测结果作为约束,并利用小波分解、纹理特征图像处理等手段,对研究区进行了地球物理场特征分析、岩石物性总结、断裂信息提取和磁性基底埋深反演计算,同时对南黄海海域三条剖面进行了广义逆重磁数据拟合反演和地震P波速度成像.根据地球物理数据的各种处理结果,提出在南黄海西部存在一条串珠状地球物理线性构造带,并将这条NNW向断续延伸的构造带称之为南黄海西缘断裂带.该断裂带延伸长,断裂两侧前新生代地层差异较大,有可能是造成陆区和南黄海南部盆地区前新生代油气远景差异的原因之一.文中得到南黄海磁性基底埋藏分布具有“中间浅四周深”的分布特征,与地震层析成像结果相吻合.磁性基底的局部起伏和局部构造的边界断层共同控制了前新生代残留盆地的格架和残余厚度分布,反映出“区域控制局部,深层约束浅层”的规律.  相似文献